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Abstract 

We explore the metastability effects across the order-disorder transition pertaining to the peak 
effect phenomenonon in critical current density ( J c ) via the first and the third harmonic ac suscep- 
tibility measurements in the weakly pinned single crystals of 2H-NbSe2- An analysis of our data 
suggests that an imprint of the limiting (spinodal) temperature above which J c is path independent 
can be conveniently located in the third harmonic data (X3u)- 

PACS numbers: 74.60. Ge, 64.70.Dv, 74.25.Dw 
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The ubiquitous occurence of an anomalous enhancement in critical current density 
(J C (H,T)) in weakly pinned superconducting samples with the increase in applied field (H) 
or temperature (T) is termed as the peak effect (PE) . There is a widespread belief that 
the PE represents a dynamical transition (e.g., from elastic flow to plastic flow) in the driven 
vortex matter jjj U, 0]. Using STM imaging measurements on a very weakly pinned single 
crystal of the well studied 2H-NbSe2 system, Troyanovsky et al p have shown that the 
collective motion of neighboring vortices gives way to positional fluctuations of individual 
vortices just above the onset temperature (T° n ) of the PE. 

As regards the nature of the phase transition across the PE region, X. S. Ling et al j]| 
have sought to provide direct structural evidence in favor of the first order nature of the 
order-disorder transition commencing at T° n in a crystal of Nb. An independent evidence 
in support of the first order nature of the phase transition in typically weakly pinned single 
crystals of 2H-NbSe2 came via bulk transport studies 0, |t| and was confirmed via local 
scanning ac micro-Hall bar microscopy from the visualization of an interface separating 
the weaker pinned (ordered) and stronger pinned (disordered) regions across the PE. The 
first order phase transition notion appears further fortified with the recent results of Xiao 



et al who explored different (partially ordered) metastable states above the usually 
ascertained T° n . Their transport measurements essentially amount to locating higher peak 
temperatures (T p ) for superheated metastable states. The highest (limiting) value of the 
T™ ax for a given field can be seen to lie close to their spinodal temperature, T s ; above it, the 
J c is ascertained to be single valued (i.e., no history effects) and it monotonically decreases 
with further increase in temperature. 

An objective of our present work is to make a proposition for the demarcation of the 
spinodal line for the order- disorder transition in the underlying (static) vortex matter using 
the contact-less ac susceptibility measurements. One advantage of the ac susceptibility 
measurements is that the vortex array is shaken by a tiny driving force, whereas in the usual 
transport experiments, to observe the PE, the vortex array has to be driven with forces 
much higher than the critical pinning force. 

We have focussed our attention on the single crystals of 2H-NbSe2- Earlier ac suscep- 



tibility studies in them [12J, |13|, have elucidated the correlation between the quenched 



random pins, effective disorder and the metastability/history effects. The PE anomaly can 
be observed to be sharper than the superconducting transition width in a nascent pinned 
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2H-NbSe2 crystal, and history effects in such a sample could be easily swamped with a tiny 
ac driving force |l5j. However, for typically weakly pinned samples, the PE phenomenon 
prior to the T C (H) comprises two distinct fracturing steps, across both of which the metasta- 
bility effects develop prominently It is further noted [15] and recalled here that with the 
progressive enhancement in quenched disorder, there appears an additional anomalous peak 
in J C (H, T) response prior to the classical peak effect located at the edge of the T C (H) line. 
The additional peak imprints as the second magnetization peak (SMP) in the isothermal 
magnetization hysteresis (M-H) loops 17]. It suffices to state here that the two anomalies 
appear distinct and different. The circumstances where they could manifest close to each 
other and/or overlap, and the dynamical behavior across them could get admixed, shall be 
dealt with elsewhere 



The ac susceptibility experiments have been performed using a home built [19( ac sus- 
ceptibility setup, in the frequency interval 21 Hz to 211 Hz, and for ac field (h ac ) am- 
plitude lying in the range of 0.5 Oe to 2.5 Oe (rms). The 2H-NbSe2 samples are the 
crystals Y' (T,W) » 7.25K) and Z (T c (0) « 6K), studied earlier by S. S. Banerjee et al 
jl2l 3, 0,13 alongwith another crystal Z' (T c (0) ~ 6K). The isofield measurements 
have been performed in three modes: zero field cooled (ZFC) warm-up, field-cooled cool 
down (FCC) and the field-cooled warm-up (FCW). 

Fig.l shows the in-phase ac susceptibility (xLCO) data in different modes in the crystal 
Y' in a field of 15 kOe. The imprint of the PE in it is self evident: x' ~ ~~ 1 + &h ac / J C {H, T), 
where a is geometry and size dependent factor j^J. A pertinent thing to note in these data 
is the distinction between the two field cooled warm-up responses. The two initial FC states 
(FC1 and FC2) were obtained while cooling down in a dc field with superimposed ac field h ac 
kept on and off, respectively. When h ac of 0.5 Oe (rms) remained present, there appeared 
little difference between the initial field cooled (FC1) state and the ZFC state, and the 
subsequent warm-up responses of both of them. However, when the h ac remained switched 
off, there is a discernible difference in the initial FC2 and the ZFC state. The difference 
during subsequent warm-up appears to survive upto (or just above) the peak temperature 
T p of the PE. These data atest to the correlation between the pinning and the exploration 



of the metastable states prior to the PE 
the disordered state by the shaking effect 



12 L and the role of an ac driving force in annealing 



2l|. 



Figs. 2(a) to 2(c) show the xLCO data in a somewhat stronger pinned crystal Z 1 (i.e., 
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where the pinning effects are larger than the crystal Y') and in which the history effects 
manifest in a robust manner such that keeping h ac on or off during cool down does not 
affect the FC state. Note first the presence of two anomalous variations in xh(T), which 
we believe represent the SMP and the PE anomalies. In the xLCO data at 2.6 kOe in 
the ZFC mode, one can also identify the two characteristic steps near T° n and T p values, 
as reported earlier by Banerjee et al At higher fields, the stepwise fracturing feature 
across the PE gradually transforms into a continuous amorphization process from T° n to 
the peak temperature T p . An interesting behavior in the data at 6 kOe and 10.5 kOe (cf. 
Figs. 2(b) and 2(c)) is that the peak temperature of the PE is different in different modes, 
the highest value being in the ZFC mode, T^ c . To emphasize, the xLCO responses in Figs. 
2(b) and 2(c) appear to mimic the scenario that emerges from the critical current I C (T) data 

Hsented for different sample histories in a crystal of 2H-NbSe2 by Xiao et al (cf. Fig. 1 of 
). The x'u>{T) response for a given H is essentially dictated by J C (T). The temperature 
above which Xu(T) (or, J C {T)) becomes independent of the thermomagnetic history of the 
specimen appears to lie even above the (highest) peak temperature, i.e., T^ c . 

The pristine enunciation of the critical state model (CSM) prescribes a specific relation- 
ship between J c and the hysteretic magnetization response of a superconductor j^j. It is 



well documented 



23, 



24 



25 



261 ] that qualitative changes can occur in this relationship when 



J C {H) does not remain single valued function of H. For instance, the minor hysteresis loops 
display anomalous characteristics, like, the asymmetric shape, excursions beyond the enve- 
lope M — H loop, etc. j3, when the J C (H) turns path dependent. In the back drop 
of these observations, it is instructive to examine the response of the third harmonic of the 
ac susceptibility across the SMP and PE regions, i.e., from below the onset temperature of 
first of two anomalous variations in J c upto the irreversibility temperature (T irr ), where the 
(bulk) J c ceases. To be specific, consider the cool down of a weakly pinned type-II super- 
conducting sample from above T c to below its T irr , where the finiteness of J c would result 
in a non-linear magnetization response which could generate a measurable third harmonic 
signal in ac susceptibility measurements. Such a third harmonic signal would be expected 
to follow the increase in J C {T) for a given H, as per prescription of the CSM for path in- 
dependent J C (H,T), as {T irr — T) increases. The onset of the history dependence in J C (H) 
could compromise the above stated notion, arising from the applicability of the CSM. We 
continue to explore below, the limit of the path dependence in J C {H) via the observation of 
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the history dependence in the first harmonic ac susceptibility data, and compare it with a 
specific feature noticeable in the temperature dependence of its third harmonic data. 

Fig. 3(a) focuses attention onto the difference plots showing (x'w FC — x'J CW ) an d (xif FC — 
x'J CC ) a t 10.5 kOe. The inset panel in Fig. 3(a) shows a portion of the data on an expanded 
scale to facilitate the marking of the limiting temperature T*, at which the larger of these 
differences vanish, and the xLCO response becomes path independent. Fig. 3(b) presents 
the plot of x'if C (T) and the difference plot, x' 3 Z f C ( T ) - x'L CW (T), in a field of 10.5 kOe. 
For the plots in Fig. 3 (b), we have identified the value of the limiting temperature T* 
determined from Fig. 3 (a) and those of T°mp^ c \ Tp and T^ c evident from Fig. 2 (c). 
Note that the x'0f C {T) data shows multiple undulations prior to reaching the limiting value 
T*. The first of these undulations appear to coincide with T°mp , where the anomalous 
variation in J c commences. We recall that the association of an enhancement in the non- 
linear response with the occurence of anomalous variation in J c is well documented in the 
literature jj. 

Next, we draw attention specifically to the behavior of x'zu j us t above T*. Warming up 
from the low temperature side, as the temperature crosses the limit T*, X^C^) is seen to 
monotonically decrease and vanish at the irreversibility temperature, T trr (< T C (H)). There 
does not appear any simple correspondence between the xL an d Xsu f° r Tsmp^^ < T <T* 
(cf. Figs. 2(c) and 3(b)). While Fig. 2(c) shows that xLCO monotonically decreases above 
Tp^ c , reflecting the collapse in J c above the peak position of PE, xk^T), on the other hand 
in Fig. 3 (b) appears to enhance between T^ c and T*. xLX^) is seen to turn around only 
above T* and follow the J C (T) thereafter (see Fig. 3(b)). As stated earlier, the signal 
can be related to the leading non-linear term in the prescription of CSM [22]. It is therefore, 
not fortuitous that an imprint of the limit of the path independence in J c is present in the 
response. 

To establish the assertion on the limit of the path independence in J c in the xkX^) 
data, we have examined the above stated behavior at different fields and in different weakly 
pinned crystals of 2H-NbSe2 and found one to one correlation between T* determined from 
x'w(T) data in different thermomagnetic histories and the limiting temperature above which 
X3u>(T) monotonically decreases. For instance, the Figs. 4(a) and 4(b) show xkX^) data 
in fields of 3.6 kOe and 5 kOe in crystals Z' and Z of 2H-NbSe2, respectively. The T° n , 
T^ c and T* values (determined from difference (x' ZFC — x'J CC ) plots) have been marked 
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for each of the curves in Figs. 4(a) and 4(b). It is evident that T* represents the limiting 
temperature above which xLjCO monotonically decreases. At T < T*, the modulations in 
X3oj(T) display complex behavior, dictated by J C {T) in different thermomagnetic histories 
and the h ac value in which the XjjwCO data are recorded. The T*, however, does not appear 
to vary in any noticeable manner with the amplitude of h ac (all data not shown here). 

It is useful to explore the correlation between the X3w(T) and the noise signal in xLCO 
jl^ l. which can be easily recorded using a Lock-in amplifier (Stanford Research Inc. Model 
SR 850) having a flat band filter option (^tJ. The above stated noise signal is believed to 
measure the fluctuations in xi){T) an d it is argued lfj to reflect the possibility of trans- 
formations amongst metastable states accessible from a given mode (ZFC/FC). Fig. 4(c) 
depicts the plot of noise in x'^ FC (T) in a field of 5 kOe and in h ac = 0.5 Oe (rms). We 
have marked the value of T* alongwith the identification of corresponding values (in h ac of 
0.5 Oe (rms) ) of T° n and T p in Fig. 4(c). It is indeed not a coincidence that the noise 
signal recedes to the background value as T — > T*. Taking cue from earlier studies of noise 

fin n 

in transport experiments |28L I23| . Banerjee et al [16] had argued that the increase in noise 
at T = T° n (= T p i in Ref. [16j) reflects the possibility of enhancement in transformations 
amongst coexisting |l0| metastable states in a fractured (partially disordered) vortex solid. 
The setting in of the decrease in the noise signal as T — > T p was considered to imply the 
effect of phase cancellation of a large number of incoherent fluctuations as the vortex matter 
moves towards the fully disordered state. In such a framework, the possibility of coexistence 
of ordered pockets embedded in the disordered medium would cease as T — > T*, and the 
noise signal would reach the background value. In the context of our present results, the 
T*(H) values represent the notion of spinodal line, T S (H) jll| . 

To summarize, we have collated together in Fig. 5 the different (reduced) field-temperature 
(h, t) values for the crystal Z'. It includes in it the t* (= T*/T c (0)) values for the other 
crystals (Y 7 and Z) of 2H-NbSe2 as well. For the sake of completeness, we have also chosen 
to depict in it the values of h pj„A t) (= H p i ateau (T) / 'ifJ^(O), where H^iO) = 44 kOe) in the 
crystal Z', as determined [r! llJ] from the normalized plots of the critical current density 
vs reduced fields (data not shown here). To recall [r3], H p i ateau represents the limiting field 
below which the collective pinning regime gives way to the small bundle pinning regime. It 
seems appropriate to identify the (h, t) space between h p i at {t) and the tf^^h) as the Bragg 
(elastic) glass region (3^]. Above the t*(h) line, where the metastability effects cease, the 
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vortex matter exists in the pinned amorphous phase. 

The parameter space in between the elastic glass and the pinned amorphous state appears 
;o belong to the co-existing weaker pinned (ordered) and stronger pinned (disordered) regions 
lol |. At this juncture we draw attention again to the xL data in Figs. 2(c) (see also, data 
in Fig. 3(a)). These data reveal that as the temperature exceeds T™ , the weaker pinned 
vortex solid in the ZFC warm-up mode starts to transform towards somewhat stronger 
pinned state (presumably with the injection and survival of disorder into the ordered state 8|, 
|9j) upto the peak temperature of the SMP anomaly. On the other hand, the transformation 
that sets in at T™ p for the stronger pinned vortex solid state in the FCW mode, takes 
this state rapidly towards the weaker pinned (ordered) state upto the onset position of the 
PE. It is indeed only above T° n that the shift towards the disordering commences for both, 
the ordered state in the ZFC mode and the disordered state in the FCW mode. While 
cooling down, the supercooling commences as temperature is lowered below T*, and it could 
appear to get further fortified as the temperature is lowered down across the SMP region 
(cf. Fig. 2(c)). We believe that in the co-existing region, somewhere the balance shifts 
from the dominance of ordered regions to that of the disordered regions, in the sense that 
transformations triggered by rise in temperature and/ or imposition of a large ac drive push 
a given state either towards a more ordered or, a more disordered state. We are tempted 
to suggest that the onset temperatures of the PE anomaly determine the said 
crossover boundary, which partitions the coexistence (h, t) space into region I 
(t°m P < t < tp n ) and region II (t° n < t < t*) (see Fig. 5). In terms of the concept of 
the stationary state [2^, T° n represents the limiting temperature above which the disordered 
regions represent the stationary state. 

Lastly, we attempt to make a contact with the theoretical results of Li and Rosen- 
stein (LR) [31] on the spinodal line, following the procedure of Ref. [ill. In terms 
of the dimensionless scaled variables, LR theory gives the spinodal line as, aLn(t,h) = 
— [(7r 2 Gz)/2]^ 1 / 3 [(l — t — h)t~ 2 / 3 h~ 2 / 3 }, where ar,Fi(t,h) = —5, and the rest of symbols have 
their usual meaning and representative values [22J, || 
take #£(0) = 44 kOe, e = 0.3, f = 8.6nm, A = 135nm 
line (solid line) satisfies our experimental data. The choice of H^(Q) = 52 kOe |5( and 
39.7 kOe, for crystals Y' and Z, make the data points for these crystals also to satisfy the 
same spinodal line 



24, 
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la]. For sample Z', if we 



361 ] . the theoretical spinodal 
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The vortex matter below the spinodal line in the phase diagram is typified by its in- 
homogenous nature of admixed weaker pinned (ordered) and stronger pinned (disordered) 
regions. The relative weights of these regions and their topographic distribution in different 
thermomagnetic histories elucidate the extent of measurable path dependence in J C (H,T). 
The presence of a driving force (ac field, transport current, electric field (^r), etc.) can 
affect the relative weights of stronger and weaker pinned regions in a dynamical manner, 
and thereby influence the observed macroscopic magnetization response between the onset 
temperature of the anomalous variation in critical current density and the limiting T* value. 
The peak temperature in anomalous variation in J c across the PE, is no long er viewed as 
the limiting temperature above which the history effects in J c cease Xiao et. al. 



11 1 have stated that the peak temperature of the PE gets determined by the competition 
between the progressive nucleation of stronger pinned regions (which has the tendancy to 
enhance the sample averaged J c value) and the decrease in pinning strength with increase in 
temperature of both the stronger and the weaker pinned regions. To this, we may add that 
the non-linear responses that can independently arise from admixed stronger and weaker 
pinned regions undergo a qualitative change above T*, where only the homogenous stronger 
pinned regions (disordered phase) survive. The T* value thus is expected to be independent 
of the amplitude of the ac driving force, whereas the peak temperature (and/or the history 
dependence in magnetization) value could depend on h ac . 

To conclude, our studies have shown that the PE is a very sharp transition in 2H- 
NbSe2 samples with very weak pinning, however, superheating/supercooling effects across 
PE are difficult to explore in such samples in presence of a driving force. In somewhat 
stronger pinned crystals of 2H-NbSe2, where the metastability effects manifest in a promi- 
nent manner, an imprint of the spinodal temperature can be conveniently located in the 
third harmonic data. 
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FIG. 1: In-phase ac susceptibility (xL) data at H = 15 fcOe(||c) and in h ac of 0.5 Oe (rms) in 
the crystal Y 1 (T c (0) « 7.25K) of 2H-NbSe2 for different thermomagnetic histories, as indicated. 
While obtaining the FC1 and FC2 states, the h ac remained on and off, respectively. 
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FIG. 2: xLO^) data (^ac = 1-5 Oe(rms)) for different thermomagnetic histories at the dc fields 
(H\\c) indicated in the crystal Z' (T c (0) « 6.0K) of 2H-NbSe2- The anomalies corresponding to 
the SMP and the PE have been identified and the respective positions of , T° n , T^ c and T c 
have been marked. 



13 



3.5 4.0 4.5 



i i i i i i i i i I i i i i i i i i i I i i i i i i i i i I i i i i i i i 




' 4.3 | 4.4 ,4.5 | }T* I 



I I I I I I I I I | I I I I I I I I I | I I I I I I I I I | I I I I I I I 




I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 
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T(K) 

FIG. 3: The panel (a) depicts the difference plots {Ax'J in the ZFC/FCC and ZFC/FCW modes 
at 10.5 kOe. The inset in panel (a) shows the data on an expanded scale to help locate the T* 
value. The values of T°™ , T° n and the T^ c identified from Fig. 2(c), have also been marked in 
the panel (a). The panel (b) shows the in-phase part of the ac susceptibility at 3lo, X3uj(T), in a 
dc field of 10.5 kOe and h ac of 2.5 Oe (rms), in the ZFC mode. In addition, it also displays the 
difference plot, A X ' 3u) {T) = X^ C {T) - xI™( T )- The value of T * determined from Fig. 3 (a) and 
those of T°mp^ c \ Tp n ^ z ^ c ^ and T^ c determined from Fig. 2(b), have been marked in the plots 
depicted in the panel (b). 
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FIG. 4: The panels (a) and (b) show X3uj(T) data at 3.6 kOe and 5 kOe for h ac = 2.5 Oe (rms) 
in the crystal Z' and Z, respectively. The panel (c) displays the noise in xL(^) recorded with a 
wide band filter at 5 kOe with h ac = 0.5 Oe(rms) and at a frequency of 211 Hz in the crystal Z. 



The positions of the respective T° n ( & T p i 
different panels. 



ld |) and T p values (see text) have been marked in the 
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t( = T/T c (0)) 



FIG. 5: A magnetic phase diagram in a typically weakly pinned single crystal Z' of 2H-NbSe2 
drawn in terms of reduced field (h) and temperature (t). The theoretical spinodal line based on 



the work of Li and Rosenstein 



3 it ] has also been drawn solid 



through various data sets are to guide the eye. Following Ref. 



ine. The rest of the lines passing 



111 ], we have also included the data 



related to the limiting (spinodal) temperature in the other 2H-NbSe2 crystals (Y' and Z). 
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